How do magnetic microwires interact magnetostatically? 
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The magnetostatic interaction between two ferromagnetic microwires is calculated as a function 
of their geometric parameters and compared with those measured through magnetic hysteresis loops 
of glass-coated amorphous Fe77.5Si7.5B15 microwires. The hysteresis loops are characterized by well- 
defined Barkhausen jumps corresponding each to the magnetization reversal of individual microwires, 
separated by horizontal plateaux. It is shown that the magnetostatic interaction between them is 
responsible for the appearance of these plateaux. Finally, using the expression for the magnetostatic 
interaction is trivial to obtain the interacting force between microwires. Our results are intended to 
provide guidelines for the use of these microwires with technological purpose such as the fabrication 
of magnetic sensors. 



I. INTRODUCTION 



In the last decades, soft magnetic materials have 
been deeply investigated. Besides the basis scientific 
interest in their magnetic properties, there are a great 
deal of technologic interest due to their use in sens- 
ing applications, particularly in the fields of automo- 
tive, mobile communication, medical and home appliance 
industries*^— Moreover, these materials are very promis- 
ing for spintronic devices in magnetic recording medial 
Two types of soft magnetic microwire families are cur- 
rently studied: in- water quenched amorphous wires with 
diameters of around 120 /jmp and quenched and drawn 
microwires with diameters ranging from around 2 to 
20 /<m^ covered by a protective insulating glassy coat. 
Bi-stable microwires are characterized by square-shaped 
hysteresis loops defined by the abrupt reversal of the 
magnetization between two stable remanent states^ 

Although an array composed of a few ferromagnetic 
wires could in principle seem a quite simple problem to 
study and model, it is striking to notice how complex 
this problem can turn out to be. Effects of interparti- 
cle interactions are in general complicated by the fact 
that the dipolar fields depend upon the magnetization 
state of each element, which in turn depends upon the 
fields due to adjacent elements. Therefore, the modelling 
of interacting arrays of wires is often subject to strong 
simplifications like, for example, modelling the wire us- 
ing a one-dimensional modified classical Ising mode h 1 ^ 11 
Zhan et ai^ used the dipole approximation including 
additionally a length correction. J. Velazquez and M. 
Vazque a 13 i 14 considered each microwire as a dipole, in a 
way that the axial field generated by a microwire is pro- 
portional to its magnetization. Nevertheless, this model 
is merely phenomenological since the comparison of ex- 
perimental results with a strictly dipolar model shows 
that the interaction in the actual case is more intense. 
They have also calculated the magnetostatic field and ex- 
panded it in multipolar terms^ showing that the non- 
dipolar contributions of the field are non negligible for 
distances considered in experiments. Besides, the mag- 
netostatic interaction energy between two magnetic ele- 
ments of arbitrary shape was derived within the frame- 
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work of a Fourier space approach by Beleggia et 
Recently, a detailed study of the magnetostatic interac- 
tion between two parallel wires placed side by side has 
been shown i 18 i 19 In spite of the extended study of the 
dipolar interactions, the effect of magnetostatic interac- 
tions on the hysteresis loops, due to the vertical displace- 
ment of the wires, has not been studied yet. 




FIG. 1. (a) Geometric parameters used for the individual wire 
description, (b) Micrograph image that shows a glass coated 
amorphous microwire. 

The purpose of this work is to develop an analytical 
model for the full long-range magnetostatic interaction 
between two microwires exploring the possibility of vary- 
ing the magnetic coupling as a function of the wires posi- 
tion. The geometry of the wires is characterized by their 
external radii R and length L, as depicted in Fig. 1. The 
separation between the wires is written in terms of the 
interaxial distance, d, and the horizontal separation, s, 
as depicted in Fig. 2. Our model goes beyond the dipole- 
dipole approximation and lead us to obtain an analytical 
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expression for the interaction in which the lengths and 
radii of the wires are taken into account. We focus on 
the interacting force in pairs of interacting wires, as a 
function of the distance between them, in order to gain 
insight on the understanding of the role of interactions 
on the plateaux that appears in the hysteresis loops. 
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FIG. 2. Relative position of interacting wires: d is the inter- 
axial distance and s is the horizontal separation. 



II. EXPERIMENTAL METHODS 

The experimental meaurements have been performed 
in glass coated amorphous bi-stable magnetic microwires 
with nominal composition Fe77.5Si7.5B15, with a satura- 
tion magnetization Mo = 1.2 x 10 6 A/m, radii of R = 6.5 
//m, and the thickness of the glass coating of t = 3.3 /jm. 
They are fabricated by means of Taylor-Ulitovsky tech- 
nique by which the molten metallic alloy and its glassy 
coating are rapidly quenched and drawn to a kind of com- 
posite microwire. Two samples with length of L = 7 mm 
each were cut from a larger wire and fixed parallel by 
means of vacuum grease on a glass holder, with a sepa- 
raton of d = 45 /xm between the wires. 

The morphology of microwires was investigated by 
a retro-optical microscope (Olimpus). The hysteresis 
curves were measured on a specially designed vibrating 
sample magnetometer (VSM), inserted within a pair of 
Helmholtz coils. These coils have sufficient field to satu- 
rate the wires and present the advantage of field homo- 
geneity and the absence of remanent fields. We will focus 
our attention on measurements performed at room tem- 
perature because a low temperature there is a change in 
the domain structure of the microwires, probably owing 
to the increasing internal stresses induced by the differ- 
ent thermal expansion coefficients of the ferromagnetic 
alloy and the covering glass J£ 



which L R so it is reasonable to assume an ax- 
ial magnetization due to shape anisotropy, defined by 
M(r) = MqZ, where z is the unit vector parallel to the 
wire axis. The magnetostatic interaction between wires 
can be calculated from2£ 

E int = fia J M 2 (r)V(7i(r)dV- , 

where M2(r) is the magnetization of the wire 2 and U\(y) 
is the magnetostatic potential of the microwire 1. The 
expression for this potential is given by 
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Note that the first term on the left-hand side of Eq. (1) 
vanishes, because the magnetization field is constant; fur- 
thermore, in the surface integral of Eq. (1) the only 
contributions arise from the upper and lower bases of the 
wire: the upper circle is located at z — L/2 and the lower 
on at z = —L/2. Due to the symmetry of the problem, we 
used the adecuate type of the cylindrical kernel for the 
integral 21 , and after few manipulations, one finds that 
the integral expression for the scalar potential is given 
by 
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where J m is a Bessel function of first kind and m order. 

Now it is possible to calculate the magnetostatic inter- 
action energy between two identical microwires using the 
magnetostatic field experienced by one of the wires due 
to the other. The final results reads 
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Equation (3) has been previously obtained for 
nanotubes32 The general expression for the interaction 
energy between wires with axial magnetization, given 
by Eq. (3), can only be solved numerically. However, 
wires that motivated this work satisfy R/L = a <C 1, in 
which case one can use that J\ (ax) ~ ax/2. With this 
approximation, Eq. (3) can be written in a very simple 
form as 



III. THEORETICAL CALCULATIONS 
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We adopt a simplified approach in which the discrete 
distribution of magnetic moments of microwire is re- 
placed by a continuous one characterized by a slowly 
varying magnetization M(r). We consider wires for 
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Finally, the magnetostatic field can be written as a 
function of the magnetostatic interaction between the 
magnetic microwires and is given by 5 
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IV. RESULTS 
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order to compensate the approximation of consider a ho- 
mogeneous magnetization in our models It is clear that 
there is a strong dependence of the magnetostatic field on 
d. Note the good agreement between experimental dat- 
apoints and analytical results for distances smaller than 
80 /*m. The coupling rapidly decreases by increasing the 
separation between the wires, and become negligible for 
interaxial distances bigger than 80 /im, at least within our 
experimental sensitivity. It can be understood if we con- 
sider that the stray field felt by one wire due to the other 
changes considerably if the wires are not completely par- 
allel to each other. Thus, when the separation between 
the wires is big enough, small deviations may reduce the 
stray field produced by the wires. 



Figure 3 shows the hysteresis loop for two microwires 
of radii of R = 6.5 /xm, length of L = 7 mm, thickness 
of the glass coating of t = 3.3 /im, and interaxial sep- 
aration d = 45 /zm at room temperature. As observed, 
the magnetization process takes place in two steps: a 
jump near 0.5 Oe, that must be ascribed to the rever- 
sion of a microwire; and an additional jump near 1.2 Oe, 
corresponding to the reversion of the other microwire. 
The flat region before the second step corresponds to the 
magnetic field felt by one wire due to the other. 
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FIG. 4. Magnetostatic field as a function of the interaxial 
separation of the microwires. The gray dots correspond to 
experimental data and the solid line represents the values cal- 
culated analytically. Parameters: R — 6.5 /J,m, L = 7 mm, 
t — 3.3 /im, s = fim and d ranging between 45 /tm and 90 
/jm. For the analytical calculations we have used an effective 
radii of R* = 4.7 ^m and M = 1.2 x 10 6 A/m. 



FIG. 3. Hysteresis loop for two parallel amorphous wires of 
Fe77.5Si7.5B15. 



A. Dependence of the interaxial separation of the 
wires on the magnetostatic field 

Theoretical and experimental magnetostatic results 
are combined in Fig. 4. Experimental data for the mag- 
netostatic field for different values of the interaxial sep- 
aration d are depicted by gray dots and the theoretical 
prediction is represented by the solid line. For the an- 
alytical calculations we have used an effective radii of 
R* = 4.7 fim., which is smaller than the real radius, in 



B. Dependence of the horizontal separation of the 
wires on the magnetostatic field 

It is interesting to analyze the behavior of the mag- 
netostatic plateaux as the interaxis distance, d, is kept 
fixed and the horizontal separation, s, is varied. Our re- 
sults are combined in Fig. 5. Experimental data for the 
magnetostatic field of the system are depicted by gray 
dots and the solid line represents the analytical calcula- 
tion. We observe a strong decrease of the magnetic field 
as the horizontal separation between wires is increased. 
Good agreement is obtained between the measured data 
and analytical calculations. 
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FIG. 5. Magnetostatic field as a function of the horizontal 
separation of the microwires. The gray dots correspond to 
experimental data and the solid line represents the values cal- 
culated analytically. Parameters: R — 6.5 um, L — 7 mm, 
t — 3.3 /im, d — 45 /im and s ranging between mm and 3.5 
mm. For the analytical calculations we have used an effective 
radii of R* = 4.7 /im and M = 1.2 x 10 6 A/m. 



C. Interacting force 

A couple of isolated magnetic microwires can attract o 
repel each other upon approach depending on their rel- 
ative orientations. For the microwires investigated the 
interacting force, F = —\7E int , has been estimated to 
be of the other of 0.017 (for d = 90 /im and s = //m) 
to 0.0685 dynes (for d = 45 /im and s = /im) for the 
case of interaxial separation in the range of parameters 
considered. For the horizontal separation, the interact- 



ing force is in the range between 1.63 x 10~ 5 (for d = 45 
/tni and s = 3.5 mm) to 0.0685 (for d — 45 /im and s = 
mm) dynes. 



V. DISCUSSION AND CONCLUSION 

In summary, we have investigated the magnetostatic 
interaction between magnetic microwires. Using a con- 
tinuous model we have obtained a simple expression to 
model the magnetostatic interaction in these particles. 
From our calculations and measurements we can conclude 
that the magnetostatic plateaux strongly depends on the 
geometry of the system. A couple of isolated magnetic 
microwires can attract o repel each other upon approach 
depending on their relative orientations. Based on mi- 
crowires investigated the interacting force, F = —VE int , 
has been estimated to be of the other of 1.63 x 10 -5 
(for d = 45 /im and s — 3.5 mm) to 0.0685 dynes (for 
d = 45 /im and s = fim) in the range of parameters 
considered. The perfomance of such experiment could 
provide the basis for testing the different theoretical mod- 
els. Our results provide guidelines for the production of 
microstructures with tailored magnetic properties. 
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